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ABSTRACT 

Context. The Galactic Center IRS 13E cluster is located ~3.2" from SgrA*. It is an extremely dense stellar association 
containing several Wolf-Rayet and O-type stars, at least four of which show a common velocity. Only half an arcsecond north 
from IRS 13E there is a complex of extremely red sources, so-called IRS 13N. Their nature is still unclear. Based on the analysis 
of their colors, there are two main possibilities: (1) dust embedded sources older than few Myr, or (2) extremely young objects 
with ages less than lMyr. 

Aims. We present the first proper motion measurements of IRS 13N members, and additionally give proper motions of four of 
IRS 13E stars resolved in the L'-band. 

Methods. The L'-band (3.8/im) observations have been carried out using the NACO adaptive optics system at the ESO VLT. 
Proper motions have been obtained by linear fitting of the stellar positions extracted by StarFinder as a function of time, 
weighted by positional uncertainties. 

Results. We show that six of seven resolved northern sources show a common proper motion, thus revealing a new comoving 
group of stars in the central half parsec of the Milky Way. The common proper motions of IRS 13E and IRS 13N clusters are 
significantly (>5<j) different. We also performed a fitting of the positional data for those stars onto Keplerian orbits, assuming 
SgrA* as the center of the orbit. Our results favor the very young stars hypothesis. 
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1. Introduction 

The central half parsec of the Milky Way is host to a 
surprisingly high number of massive young stars (see e.g. 
iPaumard et al.l 2006 Ghez et al. 2005), organized in at 



least one disk-like structure of clockwise rotating stars 



(CWS; Genzel et al. 


2003 


Levin & Beloborodov 


2003 


Paumard et al.||2006 


)• 


Paumard et al. 


(2006 


1 also propose 



the existence of a second, less populated disk of counter- 
clockwise rotating stars (CCWS). The mechanism respon- 
sible for the presence of young stars in the strong tidal 
field of the super-massive black hole (SMBH) at the po- 
sition of SgrA* is not clear. Two most prominent sce- 
narios include star formation "in-situ" (in an accretion 
disk; Levin fc Beloborodov|2003 Nayakshin et al.|2006a ), 
and the in-spiral of a massive stellar cluster formed at a 
safe distance of 5-30 pc from the Galactic Center (GC; 
Gerhard|[200T] |McMillan fe Portegies Zwart| [20(73] |Kiml 



the former scenario seems to be favored by number of au- 
thors (Nayaksh in fc Sunyaev|2005 Nayakshin et al.|2006a 



Paumard et al.||20"06 ). Also, recent results by Stolte et al 
( 2007 ) definitively rule out the possibility that the Arches 
cluster could migrate inwards and fuel the young stellar 
population in the GC. 

A special challenge and a good observational test for 
both the above mentioned scenarios is provided by the 
existence of the IRS 13 group of sources. IRS 13E (lo- 
cated ~ 3" west and ~ 1.5" south of SgrA*) is the densest 
stellar association after the stellar cusp in the immediate 
vicinity of SgrA* and contains several massive Wolf-Rayet 



(WR) and O-type stars flMaillard et al.||2004| |Paumard 
et al.|2006 Moultaka et al. 2005[ ). It is generally considered 



Vjcrna ru||zuui| |iviciv mian jroncgics Zj waru||zuuo| |nni 
et al.||2004| |Portegies Zwart et al.||2006| ). At this point 



to be associated with the mini-spiral material (Moultaka 
et al.|[2005 Paumard et al.|[2004 l and is probably bound 
since four out of seven identified stars show highly corre- 
lated velocities flMaillard et al.|2004| |Schodel et al.||2005 l 
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For both of the above mentioned scenarios of star for- 
mation at the GC there are several issues when dealing 
with IRS 13E. In principle, such a cluster could have been 



formed in an accretion disk (Milosavljevic & Loeb 2004 



Nayakshin & Cuadra 2005 1. However, in numerical simula- 



tions of star forming disks, fragmentation of a disk cannot 
produce such a dominant feature (Nayakshin et al. 2007). 
In light of the cluster in-fall scenario, an intermediate- 
mass black hole (IMBH) was proposed to reside in the 
center of the cluster ( Maillard et al.|[2004 1 . The existence 
of an IMBH makes the process of cluster in-spiral much 
more efficient in terms of increasing dynamical friction 
and thus allowing the most massive stars that reside in 
the center of a very massive (>10 6 M Q ) stellar cluster to 
reach the central parsec of the Galaxy within their life- 



times (Hansen & Milosavljevic 2003| Berukoff & Hansen 
[2006] [Portegies Zwart et al.||2006| but see |Kim et al.||2004 



for a characterization of the problems with this hypoth- 
esis). However, the presence of the IMBH in IRS 13E is 



disputed. Schodel et al. (2005 1 analyze the velocity disper- 



sion of cluster stars and derive that the mass of such an 
object should be at least 7000 M Q , making its existence 
implausible. Both the X-ray (Baganoff et al. 20031 and 



radio (Zhao & Goss 19981 source at the position of IRS 
13E can be explained by colliding winds of high-mass los- 
ing stars ( |Coker et al!|2002| |Zhao fc Gos"s||1998[), without 



the need for any unusual object. Paumard et al. (20061 
suggest that, in the case that IRS 13E is associated with 
the Bar of the mini-spiral, the mass of the stellar content 
would be high enough to keep the cluster from disruption. 

A second comoving group of stars, called IRS 16SW 
and located 1.9" in projection from SgrA*, was reported 
by Lu et al. ( 2005 1 . Four of the five members of the 



IRS 16SW comoving group are spectroscopically identi- 
fied as young massive stars, indicating ages comparable 
to those of the IRS 13E cluster. However, unlike IRS 13E, 
IRS 16SW shows only a slight stellar number density en- 
hancement. 

Approximately 0.5" north of IRS 13E, a small clus- 
ter of unusually red compact sources has been reported 
(IRS 13N; |Eckart et al.|2004[ ). Eight sources have been re- 
solved and labeled a through -q as shown in Figs. 1 and 2 



of Eckart et al. ( 2004 1 . A strong infrared excess is due to 
the emission of warm (T -1000K) dust (|Moultaka et~aT 



20051. The authors propose two possible explanations for 
the nature of IRS 13N: (1) objects older than few Myr 
and similar to bow-shock sources reported by |Tanner et al.| 
( 2005[ ) or (2) extremely young stars (0.1 - 1 Myr old). The 
latter scenario implies more recent star formation than 
what has been assumed so far for the GC environment. 

Here we present first proper motion measurements of 
the IRS 13N sources and additionally present proper mo- 
tions of four of IRS 13E sources, including one for which no 
proper motion has previously been published. All proper 



Table 1. Proper motions of IRS 13N sources observed in 
L'-band. 



Q 

P 

7 
6 
e 

c 

13!Sj] 



Aq 
^2/70 
-2.89 
-3.10 
-2.92 
-2.88 
-3.16 
-3.11 
-3.01 



~Ksf 

-1.48 
-1.25 
-0.99 
-0.85 
-1.02 
-0.82 
-0.66 
-0.93 



VR.A. 

44 ± 8 
- 6 ±20 

-114 ±14 
-33 ±13 
-58 ±15 

-134 ± 9 
-91 ±22 
-73 ± 7 



VDec. g 

147 ±13 
228 ±15 
298 ± 7 
248 ± 9 
306 ±13 
333 ± 6 
323 ±39 
289 ± 8 



"relative to SgrA*, in arcseconds 
'all velocities are in kms -1 ; the uncertainties 
represent the la uncertainty of the linear fit 

c average proper motion of six comoving sources f3-r] 



Table 2. Proper motions of IRS 13E sources observed in 
L'-band. 



name 



Aa 



El 
E2 
E3 
E5 
13E 
13E 



-2.95 
-3.16 
-3.17 
-3.39 
-3.17 
-3.24 



VR.A. 



VDe 



-1.64 
-1.73 
-1.53 
-1.53 
-1.61 
-1.60 



-107 ±11 
-228 ± 8 
-158 ± 4 
-175 ±45 
-167 ±12 
-187 ±15 



-118 ± 6 
53 ± 9 
71 ± 3 
140 ±72 
37 ±18 
88 ±24 



"relative to SgrA*, in arcseconds 
'all velocities are in kms -1 ; the uncertainties 
represent the la uncertainty of the linear fit 

c average proper motion of four IRS 13E sources 
d average proper motion of IRS 13E without El 



motions were derived using L'-band NACO/VLT images 



2. Observations and Data Reduction 

The L' (3.8 /im) images were taken with the 
NAOS/CONICA adaptive optics assisted im- 



ager/spectrometer (Lenzen et al. 




1998 


Rousset et al. 


1998 Brandner et al. 




2002) at 


the UT4 (Yepun) at 



the ESO VLT. The data set includes images from 6 
epochs (2002.66, 2003.36, 2004.32, 2005.36, 2006.41 and 
2007.39) with a resolution of ~100mas and a pixel scale 
of 27 mas/pixel. Data reduction (bad pixel correction, sky 
subtraction, flat field correction) and formation of final 
mosaics was performed using the DPUSER software for 



astronomical image analysis (T. Ott; see also Eckart & 
Duhoux|1990t . All the images were deconvolved using the 



1 Based on observations collected at the European Southern 
Observatory, Chile 




linear Wiener filter. The absolute positions of sources in 
our AO images were derived by comparison to the VLA 
positions of IRS 10EE, 28, 9, 12N, 17, 7 and 15NE as 



given by Reid et al. (2003) 



Stellar positions were extracted with StarFindcr 



(Diolaiti et al. 20001 and transformed into the common 



coordinate system with the aid of 19 reference stars. The 
positions of the reference stars were corrected for the stel- 
lar proper motions as derived from the Kg-band images 



(see |Muzic et al.||2007[ ). In |Muzic et al.| ( |2007[ ) we have 
shown that the lower resolution L'-band data can be reli- 
ably used to obtain proper motions. 

Proper motions were derived by linear fitting of po- 
sitions as a function of time, weighted by the positional 
uncertainties. Both the error of the transformation to the 
reference frame and the error of the stellar position fitting 
contribute to the uncertainties. 

We assume the distance to the GC to be 7.6 kpc 
( Eisenhauer et al.||2005| . 



3. Results 

3.1. Proper motions 

In tables [T] and [2] we list proper motions of IRS 13N 
and IRS 13E sources identified in L'-band, respectively. In 
Fig. [T] we show proper motions of all stars superposed on 
the L'-band image. Stars (3 through n show similar proper 
motions, revealing a new comoving group of sources at the 



GC. As shown in Eckart et al. (2004), only rj can also be 
detected in the K-band. Paumard et al. ( |2006[ ) give the full 
velocity information for this star: vr.a. = (-52±28)kms~ 1 , 
v Dec = (257±28)kms~ 1 and v r = (40±40)kms~ 1 . The au- 
thors identify it as a member of a possibly existing second 
stellar disk (CCWS). 



Concerning IRS 13E, the obtained proper motion for 
El shows the biggest discrepancy with respect to previ- 
ously published data ( Paumard et al.|[2006 Schodel et al. 
2005). It also deviates significantly from the proper mo- 



tions of the other three stars. It was already noted that 
El shows the largest deviation from the systemic proper 
motion of the cluster (jSchodcl et al. 2005). Paumard et al. 
(2006) argue that this star is possibly not bound to the 



cluster. This is supported by our result. 

It is important to note that, in addition to a very dense 
stellar population in the region, a large amount of diffuse 
dust emission detected at 3.8 /xm makes the astrometry 
more challenging to perform than at shorter wavelengths. 
For this reason, StarFinder was not able to resolve all 
the previously reported K-band sources in the region. E3 
is known to be at least a double source, while here we 
present it as a single one. E4 is located close to E3 but 
is much fainter in the L'-band. Therefore it could not be 
clearly identified by StarFinder. However, we were able 
to resolve the faint source E5 in all epochs and present 
the first proper motion measurement of this star. E5 was 



proposed to be a dusty WR star (Maillard et al. 2004), 



and it seems to exhibit a proper motion consistent with 
those of E2 and E3. This implies that it may be a cluster 
member. 

From Tables [T] and [2] one can see that the average 
proper motions of the two clusters are significantly (> 5 a) 
different. Thus, it seems that two clusters do not belong 
to the same system. 

3.2. Keplerian orbit fitting 

In order to analyze our measurements we attempted to 
fit the positional data of IRS 13N stars showing similar 
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(a) 




4 2 -2-4 6 4 2 -2 -4 -6 8 6 4 2 -2 -4 -6 

R.A. [arcsec] R.A. [arcsec] RA [arcsec] 

Fig. 2. (a) the best fit orbital solutions for five IRS 13N stars (7-77); (6) orbital solutions for a single plane (i—2A°, 
f2=180°), for stars 7 - 77 ; (c) chosen la orbits for 77, with the best fit orbit colored black; NOTES: part of the orbit 
in front of the plane of the sky is colored red; crosses mark present positions of the stars; axes show the offset from 
SgrA*. 



proper motions (stars (3 to 77) to Keplcrian orbits, assum- 
ing that the gravitational potential is dominated by the 
SMBH at the position of SgrA*. For details of the fitting 
procedure and definition of parameters see e.g. |Eisenhauer| 
et al.l (20031. We assumed the distance and the mass of 



SgrA* to be 7.6 kpc and 3.6 xlO 6 M , respectively. The ve- 



locity of SgrA* was assumed to be zero ( Reid & Brunthaler 



20041. Under the assumption that the stars are really on 



bound orbits around SgrA*, it is clear that our data cover 
only a small part of such an orbit. We also lack the radial 
velocity information for all of the stars. As a consequence, 
the x 2 minimization gives a very wide span of orbital pa- 
rameters for which x 2 <x 2 m m+l, even when fixing most 
of the free parameters (the mass of and the distance to 
SgrA*, as well as the position and the velocity of the cen- 
ter of the orbit). Therefore we repeated the orbital fitting 
for each star introducing the radial velocity data as given 
in iPaumard et al. (20061 for 77. We varied the inclination 



angle i between 0° and 90° in steps of 1°, f2 (position angle 
of the ascending node) between 0° and 180° in steps of 20° 
and uj (longitude of periastron) between 0° and 360° with 
a step size of 20°. The results are given in Table [3] We list 
the best fit parameters i and f2, and give the la range for 
the inclination. We also list the range in eccentricities (e) 
and semi-major axes (a) resulting from the fits for which 
there is some combination of angles (z, f2, ui) that gives 
a fit better than la For all the stars except £, the 
best fit reduced \ 2 value ranges between 0.6 and 1.5. For 
£ the best fit is obtained with % 2 =3.9, probably implying 
that either positional uncertainties are underestimated, or 
that assuming the radial velocity of 77 is a poor assump- 
tion. In the last three lines of Table [3] we give the best-fit 
parameters for the Keplerian orbits fitted onto the aver- 
age positions of the stars within the IRS13N and IRS13E 
comoving groups. For IRS 13E we perform the fitting pro- 



2 a fit with x 2 Ky? min + 1, for a corresponding star. For a 
convenience, we refer to all orbits that satisfy this condition as 
la orbits. 



cedure twice: first including the star El and afterwards 
considering it a non- member. 



Paumard et al. (20061 conclude that IRS 13E is on a 



fairly eccentric orbit, with e>0.5, which we confirm. For 
all IRS 13N stars, with the exception of /3, the range in 
eccentricities of la orbits is restricted to fairly low values 
(< 0.4), in agreement with the |Paumard et aL] ( |2006[ ) value 
listed for 7;. Note that the corresponding semi-major axes 
also result fall within a restricted range of values. In con- 
trast, both parameters (e and a) are poorly constrained 
for the star (3. The size of the best fit orbit of (3 is an or- 
der of magnitude higher than those of the other five stars. 
Since the orbit of (3 is radically different from orbits of all 
the other stars, we note that it probably should not be 
regarded as part of the co-moving group. Therefore we do 
not show its orbit in Fig. [2] Since 77 is the only star with 
known radial velocity, we took its orbital parameters and 
attempted to fit other stars onto the same orbit. We note 
that it is not possible to fit all other stars onto exactly the 
same orbit (same i, £1 and u>). Therefore we fix i and £1 
and let uj vary. This confines the stars to the same plane, 
within which they can still have different orbits. For three 
stars (7, e, Q there exist la orbits with given i and f2, 
the star S can be reasonably (but still not within la) fit- 
ted, and (3 fails completely to be fitted onto this plane 
(Fig. [2^). To get a feeling about the uncertainties of our 
analysis, we plot several la orbits for the same star (77) in 
Fig#. 

The implications of these results are further discussed 
in Section IX2I 



4. Discussion 

4.1. 3-dimensional position 



Paumard et al. (20041 report the morphology and ra- 



dial velocities of the mini-spiral gas based on He I and 
Br 7 observations. At the position of the IRS 13 sources, 
the Northern Arm and the Bar of the mini-spiral are 
overlapping along the line of sight, with the Bar being 
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Table 3. Results of the Keplerian orbit fitting 



5 



P 
7 
<5 
e 

c 

n 

IRS 13N 
IRS 13E 
IRS 13Er¥ 



79 
19 
48 
33 
18 
24 
28 
83 
70 



180 
180 
180 
180 
180 
180 
180 
100 
120 



53 
4 - 
39 
17 

- 

1 - 
14 
80 

14 - 42: 



80 
67 
71 

56 
30 
45 

49 
■ 85 
58 - 



80 



- 
0.17 
0.15 
0.1 
0.1 
0.1 
0.1 
0.5 
0.5 



"best fit parameters (x 2 ~ X 2 min) 

6 all inclination values for which x 2 <X 2 min + l 



0.9 
0.21 

- 0.4 
0.3 
0.3 
0.3 
0.2 
0.9 
0.8 



4.0 
2.6 

2.7- 
2.5 
3.2 
2.8 
2.8 
12.5 - 

2.3 - 



Calculated for (i, fi, 
d average orbit 
including El 
hiot including El 



combinations for which \ <x min+1 



180.0 
4.0 
15.0 
7.2 
4.7 
5.7 
5.7 
145.0 
66.0 



placed further away from the observer. The authors sug- 
gest that IRS 13E must be either embedded in the Bar, or 
very close to it. Also, the radial velocities of the IRS 13E 
stars are similar to the radial velocities of the Bar mate- 
rial at this position ( Paumard et al.|[2004 1 . Spectroscopic 



observations of the entire IRS 13 region (Moultaka et al. 
20051 indicate a close spatial correlation between the stel- 



lar sources and the surrounding material. Higher water 
ice and hydrocarbon absorptions in the northern part of 
IRS 13, as well as a redder continuum emission, suggest 
that the IRS 13N sources are more embedded in the ISM 
than the rest of the complex. The radial velocity of ry is 
consistent with the radial veloci ty of the Bar at this po- 
sition (between and 50kms _1 , Paumard et al.|[2004 ). If 
the northern sources that show similar proper motions also 
have correlated radial velocities, it would be reasonable to 
suggest that the entire complex is indeed associated with 
the Bar. According to our orbital fits, IRS 13N is located 
a few equivalent arcseconds behind SgrA*, which nicely 



agree with the view of Paumard et al. ( 2004 1 about the 



Bar being located somewhat behind SgrA*. This would 
imply that the IRS 13E and IRS 13N stellar complexes are 
spatially close. The existence of two comoving clusters, 
each moving in a completely different direction, poses a 
lot of challenges for star formation scenarios at the GC. 



4.2. Is IRS13N a bound system? 

As can be seen from the orbital analysis, the best fit orbits 
for IRS 13N stars have rather different orbital parameters. 
Observing six stars on different orbits exactly at the mo- 
ment when they appear to be very close in projection and 
to have similar proper motions is not very probable. There 
is no single orbit onto which all of the stars would fit with 
a good x 2 , but at least four of the stars can be confined to 
the same plane. In this case the stars are again on differ- 
ent orbits and consequently have different orbital periods, 



spanning the range from 1000 to 3000 years. This again 
implies that the present arrangement is temporary. Also, 
it is curious that the moment at which we observe them co- 
incides exactly with the presence of dense gas at their po- 
sition. Interestingly, their most likely common plane seems 
to coincide with the plane of the counter-clockwise mov- 
ing stars (CCWS; [Paumard et aL|2006[ ). The existence of 
the CCWS disk is still a matter of debate. According to 



Paumard et al. (20061, the disk is sparsely populated (12 



stars), with rj belonging to it. The assumption that all 
IRS 13N stars belong to the CCWS significantly increases 
the population of the disk, and at the same time weakens 
the claim that the CCWS is essentially in non-circular mo- 



tion. In the analysis of Paumard et al. ( 2006 1, seven out of 



twelve CCWS stars are on eccentric orbits (e>0.4), three 
of those belonging to IRS 13E system. We find that five of 
the IRS 13N stars are on low-eccentricity orbits. 

The analysis of la orbits for one star (Fig. |2p) shows 
that the uncertainty of semi-major axes (o), as well as 
that of orbital periods (P), is of the same order as the dif- 
ferences between a and P of each star. Therefore, at this 
point we still cannot exclude that IRS 13N stars are indeed 
orbiting together. In the following we consider the possi- 
bility that the cluster is bound. We calculate the velocity 
dispersion of stars as 



JV 



i=l 



' error 2 {v x .i) + error 2 (v y>i )]/[2(N - 1)] 



where a pm is the dispersion of the measured proper mo- 
tions and the second term removes the influence of the 
proper motion measurement uncertainties. 

The mass estimate from the velocity dispersion of 
stars gives ~3300M Q , an unrealistically high mass for 
such a cluster. For a comparison, we note that |Paumard| 
et al. (20061 estimate the total stellar mass of IRS 13E 
to be ~1OOOM0. We note that even if the uncertainties 
in proper motion measurements were underestimated, the 
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final mass of the cluster calculated from the velocity dis- 
persion would remain high. For instance, a hypothetical 
underestimation of our uncertainties of 30% would result 
in ^8% lower velocity dispersion and the mass of 2700 M Q . 
The Hill radius 



. / to \ s 



for a cluster of m = 3300M Q , with a =(2.8 -5.7)" and 
e = 0.1 -0.3 is rnui ~ 0.17" - 0.34" , with the average value 
remarkably close to the observed cluster radius of ~ 0.25" . 
The other way around, if we set rum ^ 0-25" , the maxi- 
mum a and minimum e, the minimum mass needed to 
keep the cluster from disruption would be 1250M Q . 
The relaxation time for a system of N stars with the mean 
mass to* is given by 



N 



X t r 



where t cross =R/v and A=i?v 2 /Gm*«iV (Binney & 
Tremaine|1987[ ). Also, it is important to note that tidal in- 
teractions tend to shorten relaxation timescales. A cluster 
with the IRS 13N velocity dispersion would have a t re i ax 
of the order of the IRS 13N orbital period only if it was 
consisted of 500-1000 stars. For m^3300M Q contained in 
500-1000 main sequence stars we get that stars should be 
as bright as tox=17-18 (assuming the extinction Ajf «3, 
Schodel et al. 20071. These stars would be observed. For 
a lower mass limit of 125OM0, we get to^=19-21, exactly 
at the limit of NACO K-band data. Therefore it seems 
that the system could hardly be bound and survive for a 
significant time in the orbit around SgrA*. 

If IRS 13N stars are not bound, the observed veloc- 
ity dispersion probably indicates that the cluster is in the 
process of dissolution. The velocity dispersion will sub- 
sequently diminish the stellar surface number density of 
the cluster, which would within only a few hundred years 



reach the background value given by Schodel et al. (2007) 



4.3. The nature oflRS13N 



account that the Galactic Center is an extremely com- 
plex environment and that any star formation process that 
takes place there is completely different from the "nor- 
mal" star formation via the fragmentation of a molecu- 
lar cloud. Timescales for star formation could be signifi- 
cantly shorter. The free-fall time in a collapsing cloud at 
the distance of the young stars in the central half parsec 
around SgrA* is of the order of 100-1500 yr. Also, simu- 
lations of star formation in a gaseous disk around Sgr A* 



by Nayakshin et al. ( 2007 ) indicate that the timescales for 



star formation are fairly short, of the order of several thou- 
sand years. This is comparable to the orbital timescale of 
IRS 13N. 

A serious obstacle for forming stars at the GC is that 
the gas densities in the central parsec are way too low for 
the gas to be able to form stars. Even the highest density 



estimates for the circumnuclear disk (CND; Christopher 
et al.|2005 1 are several orders of magnitude lower than re- 
quired to fulfill the criterion for Jeans instability. However, 
with the aid of shocks and collisions, a small clump of gas 
in-falling to the center could be highly compressed and 
star formation within it triggered. As a matter of fact, 
the mini-spiral is a short-lived feature (tdyn^W 4 yr) and 
is in-falling right now. Furthermore, the material of the 
mini-spiral is susceptible to shocks, as discussed by e.g. 



Muzic et al. (20071. A possibility of forming stars in small 



associations has never been discussed so far for this re- 
gion. Star formation in small groups definitely would not 
be able to account for the young co-eval stellar population 
in the central parsec:, but could, in general, be possible and 
result in small associations like IRS 13N. 

However, it is clear that we still must be very careful 
when drawing a conclusion that IRS 13N cluster contains 
very young stars. In case that the cluster is indeed a part 
of the CCWS disk and that the stars were formed simul- 
taneously with other CCWS stars, then they must be sev- 
eral Myr old. The fact that they are apparently strongly 
extincted may then indicate that they are similar to the 



mini-spiral bow-shock sources ( Tanner et al. 2005 1 , but 



many times less massive and luminous (see discussion in 
Eckart et al. 2004). The question remains: why are they 



The detailed discussion about the nature of IRS 13N clustered? 



sources is presented in Eckart et al. (2004). Among other 



possibilities, the authors argue that infrared excess sources 
IRS 13N could represent objects that have colors and 
luminosities consistent with YSOs and Herbig Ae/Be 
stars. Dusty envelopes that usually surround those objects 
would then give rise to the observed strong infrared excess. 
It is striking that six of the sources that are very close 
in projection, also show very similar proper motion val- 
ues. In the previous section we argue that IRS 13N could 
hardly survive in the present arrangement for a significant 
amount of time. The indication of a dynamically young 



stellar system concurs with the Eckart et al. (2004) hy- 



pothesis of IRS 13N stars being extremely young objects. 
But why then does this object exsist at all? Timescales for 
its disruption seem to be much shorter than any timescales 
within which stars normally form. Here we must take into 



5. Conclusions 

In this paper we present the proper motion measurements 
of IRS 13N, a newly discovered co-moving group of red 
sources at the GC. We discuss the boundness of a cluster 
and a possibility that these stars are the youngest stellar 
object ever observed in this environment. Ongoing star 
formation in the GC would have significant implications 
on our understanding of this particular environment, and 
star formation in galactic nuclei in general. It has been 
assumed that the star formation in the GC occurred in 
two star bursts, ~100Myr and ^7 Myr ago (Krabbe et al. 



1995 ). The existence of objects younger than 1 Myr would 



imply that in-situ star formation in the GC and near mas- 
sive black holes in general is possible. Furthermore, this 
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would also mean that star formation in the GC must not 
necessarily be related to star bursts and that, in fact, a 
continuous star formation could go on. 

The question of the real nature of the IRS 13N cluster 
remains open. More imaging and especially high resolution 
spectroscopic data that could reveal features attributed to 
YSOs will be necessary to confirm our hypothesis that the 
stars are extremely young. 
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